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Abstract 

We have computed the fuU one-loop corrections (electroweak as weU as QCD) to processes 
contributing to the rehc density of dark matter in the MSSM where the LSP is the lightest 
neutralino. We cover scenarios where the most important channels are those with gauge boson 
pair production. This includes the case of a bino with some wino admixture, a higgsino and 
a wino. In this paper we specialise to the case of light dark matter much below the TeV 
scale. The corrections can have a non-negligible impact on the predictions and should be taken 
into account in view of the present and forthcoming increasing precision on the relic density 
measurements. Our calculations are made with the help of SloopS, an automatic tool for the 
calculation of one-loop processes in the MSSM. The renormalisation scheme dependence of the 
results as concerns tan f3 is studied. 



1 Introduction 



Existence of nonbaryonic dark matter is established through precise determination of the mean 
densities of matter in the Universe. However one does not know what constitutes dark matter even 
if the measurement of the rehc density of cold dark matter is now very precise[T]. On the other 
hand the latest observations dedicated to the detection of dark matter have recently received a 
lot of attention: the PAMELA collaboration has reported a 100 GeV excess on the ratio of fluxes 
of cosmic ray positrons to electrons^ and the ATIC balloon experiment claims a cut-off in the 
positron flux near 500 GeV 3 . The FERMIj4j and HESS^Sj observations do not confirm this data 
but still point to a deviation from the power-law spectrum. Many explanations were advocated to 
account for these results and the most optimistic and exciting one is that it could be a signal of an- 
nihilating dark matter. Supersymmetry can provide, among many other advantages, a dark mater 
candidate through the lightest supersymmetric particle (LSP) which is a neutralino, Xii if R-parity 
is conserved. Meanwhile, the search for dark matter will soon also take place within colliders, in 
particular the LHC. If dark matter is discovered among the other new particles that form a New 
Physics model, one will be able to probe its properties. One could then predict the relic abundance 
of the Universe and would constrain cosmology with the help of precision data[51[71[5 provided by 
WMAP[lJ and PLANCK^. The present WMAP accuracy on the relic density is about 10% and 
with the PLANCK satellite that has been launched recently it will reach about 2% precision. So- 
phisticated codes exist [TOl [Til HI] for the calculation of the relic density in supersymmetry with the 
inclusion of some higher order effects, essentially through some running couplings/masses or some 
effective couplings (particularly corrections to the Higgs couplings that can drastically change the 
results in the so called Higgs funnel [IIJIITS] for example). However these codes are essentially based 
on tree-level cross sections. To match the experimental accuracy, on the theoretical side we have to 
provide more precise calculations. Therefore we need to evaluate annihilation and co-annihilation 
cross sections at least at next-to-leading order. Considering the very large number of processes 
required for the evaluation of the relic density and the number of diagrams that each process in- 
volves, especially at one loop, automation of the relic density calculations especially in the MSSM 
in unavoidable. The purpose of this paper is to present some results on the one-loop calculation 
of the relic density of the LSP in the MSSM, where the dominant annihilation and co-annihilation 
channels are dominated by annihilations into gauge bosons. Beside the physics motivation for 
such scenarios, calculations of these processes involving gauge bosons are challenging. Attempts 
to include some effects through effective couplings are tricky because one has to be careful about 
maintaining gauge invariance and unitarity. A preliminary study of such scenarios has been made 
by some of usf6lll4j. In this paper we consider the case of a relatively light neutralino. Very heavy 
neutralinos with TeV masses and above will be studied in a subsequent paper especially since they 
show new interesting effects. 

The relic abundance will be derived from the assumption that it is thermally produced in the early 
stages of the universe|15|. so in a first approximation the relic density is inversely proportional to the 
thermal average cross sections, {(tv). The computations that we present here are performed with the 
help of SloopS[16, ijj J7l[T8], which is a fully automated code for the NLO calculation of any cross 
section or decay in the MSSM. Although our main interest concerns the channels with gauge bosons 
in the final state, we will find that there are non negligible co-annihilations channels with quarks 
in the final state. We calculate both the electroweak and QCD corrections to these contributions. 
There can be a host of processes contributing to the relic density for a particular scenario. In this 
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paper we calculate however the radiative corrections only to those processes which, at tree-level, 
contribute more than 5% to the relic density. We study here three different scenarios corresponding 
to three different compositions of our lightest neutralino: i) a bino like neutralino with some wino 
admixture , ii) a higgsino like neutralino, iii) a light wino like neutralino. We will also study the 
impact of different choices of the renormalisation scheme for tan/3 on the corrections. To derive 
the corrected relic density we feed our velocity dependent cross sections into micrOMEGAsfT^ for 
performing the thermal average and solving the Boltzmann equations. We will always show the 
cross sections (at tree-level and at one-loop) in terms of the relative velocity and as a guide we 
extract the s-wave and p-wave coefficients and the corrections they receive. This helps also extract 
the, one-loop, Sommerfeld[Tn] factor for some of the cross sections. In the processes we study here 
these are of QED origin. Once these one-loop Sommerfeld QED corrections are extracted we first 
subtract them before performing the all order resummation and deriving the relic density. 
The paper is organised as follows. In the next section we give an overview of SloopS and on how we 
perform the one-loop calculations, in particular summarising our renormalisation procedure. Checks 
on the calculations as concerns ultra-violet finiteness, infra-red finitcncss and gauge parameter 
independence are spelled out. The interface between SloopS and micrOMEGAs will also be presented. 
Most of the scenarios that we will be studying involve co-annihilation, we will define the effective 
cross section that includes the statistical weight. At the end of this section we present how our 
models have been defined. Section 3 considers the case of a light (about 100 GeV) mixed bino-wino 
LSP. The case of a dominantly higgsino LSP with mass of about 200 GeV is studied in Section 4. 
Section 5 covers the case of a light wino of about the same mass. The last section contains our 
conclusions and prospects. 

2 Overview of the calculation 

2.1 Set up of the automatic calculation: SloopS 

One-loop processes calculated via the diagrammatic Feynman approach involve a huge number of 
diagrams even for 2 — > 2 reactions, especially in a theory like supersymmetry. Doing full calculations 
by hand without automation is practically untractable. There exists already efficient automatic 
codes for one-loop calculations [20 l [2T | I22j. SloopS is an automated code for one-loop calculations 
in supersymmetry. It is a combination of LanHEP [23] . the bundle FeynArts|24j. FormCalc[22j and an 
adapted version of LoopToolspSl ^161 (that we will call the FFL bundle from now on). LanHEP deals 
with one of the main difficulties that has to be tackled for the automation of the implementation of 
the model file, since this requires that one enters the thousands of vertices that define the Feynman 
rules. On the theory side a proper renormalisation scheme needs to be set up, which then means 
extending many of these rules to include counter-terms. This part is done through LanHEP which 
allows to shift fields and parameters and thus generates counterterms most efficiently. The ghost 
Lagrangian is derived directly from the BRST transformations. The loop libraries used in SloopS 
are based on LoopTools with the addition of quite a few routines in particular those for dealing 
with small Gram determinants that appear in our case at small relative velocities of the annihilating 
dark matter, and even more so of relevance for indirect detection[16j. 
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2.2 Non-linear gauge fixing 

We use a generalised non-linear gauge [26t [16] adapted to the minimal supersymmetric model. The 
gauge fixing writes 

Lgf = - ^\{d^, ~ ieaA^ - lgcwPZ^)W^'+ + iS,w^{v ^Ih"^ + CjH"^ + ikG"^ + ipA'^)G+ 

- T^id.Z^^ + + ^ ~^)c0)2 „ J_(^d,A^f . (1) 

Unlike the other parts of the model Cgf is written in terms of renormalised fields and parameters. 

are the Goldstone fields. We always work with £,a,z,w = 1 so as to deal with the minimal 
set of loop tensor integrals. This implementation of the gauge fixing is very useful to check gauge 
independence for processes with gauge boson production. More details are given in|17]. 

2.3 Renormalisation 

In our code we have renormalised and implemented each sector of the MSSM. This is explained in 
details in[Sl [T71 HH]. Here we only briefiy sketch the renormalisation procedure. We have worked, 
as far as possible, within an on-shell scheme generalising what is done for the electroweak standard 
model ;2D]. 

i) The Standard Model parameters : the fermion masses as well as the mass of the W and Z are 
taken as input physical parameters. The electric charge is defined in the Thomson limit, see for 
example [20]. The light quarks (effective) masses are chosen such as to reproduce the SM value of 
a^^(M|) = 127.77. This should be kept in mind since one would be tempted to use a DR scheme 
for a, defined as M^, to take into account the fact that dark matter is annihilating at roughly the 
electroweak scale, so that a(M|) is a more appropriate choice. We should keep in mind that doing 
so would amount to correcting the tree-level cross section by about 13% for 2 ^ 2 processes. As 
we will see this running does not, most of the time, take into account the bulk of the radiative 
corrections that we report here. 

ii) The Higgs sector : We take the pseudoscalar Higgs mass Ma as an input parameter and re- 
quire vanishing tadpoles, tan/3 is defined through several schemes whose impact on the radiative 
corrections we will study: 

- a DR definition where the tan (3 counter-term is defined as a pure divergence leaving out all finite 

parts. 

- a process-dependent definition of this counter-term by extracting it from the decay — > t'^t~ 

that we will refer to as At-,- for short. This definition is a good choice for the gauge indepen- 
dence of the processes. 

- an on-shell definition with the help of the mass of the heavy CP Higgs taken as input 

parameter called the MH scheme from now on. We have reported elsewhere that this scheme 
usually introduces large radiative corrections. 

These schemes are thoroughly discussed in [17,, in particular the question of gauge invariance of 
these schemes is addressed. 
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iii) The sfermion sector : For the slepton sector we use as input parameters masses of the two 
charged sleptons which in the case of no-mixing define the R-slepton soft breaking mass, Mg^ and 
the SU{2) mass, Mg^^, giving a correction to the sneutrino mass at one-loop. In the squark sec- 
tor each generation needs three physical masses to constrain the breaking parameter Mq^ for the 
SU{2) part, Ms„, Mj^ for the R-part. See[18l for details. 

iv) The chargino/neutralino sector : For this sector we implement an on-shell scheme by taking 
as input three masses in order to reconstruct the underlying parameters Mi, M2, 11. In SloopS [IB] 
the default scheme is to choose two charginos masses and 'm^± as input to define M2 and fi 
and one neutralino mass, rn^o, to fix Mi. The masses of the remaining three neutralinos receive 
one-loop quantum corrections. In this scheme, these counterterms are |18) 

1 / Sm-+ 6in^+ 

SM2 = -779 7 (M2m2+ - fidctX) ^ + (AfafTi + - ^^detX) ^ 

Mi - \ '^i rn^+ X2 m-+ 

\ ^1 A2 

- M2,(M2 + MS2/3)^^ - MM^S2/3C2/3^ 

I ( 6m^+ Sm-+ 
= —2 rF2 - M2detX) ^ + {fiml+ - A'hdetX) ^ 

- M^ifi + M2S2/3)^3^ - M2M^S20C2p^) , (2) 

Mw */3 J 

^Mi = -^{5^x1 - Nli5M2 + '2Nl^Nl^5^l 

- 2iV*i7V*3(5Yi3 - 2Nl2Nl^5Y2^ ~ 2NliNl^5Yii - 2Nl2Nl^5Y2^) , (3) 

with detX = M2IJL — M^S2/3, Y is the neutralino mixing matrix and N its diagonalising unitary 
matrix, see |18j . Looking at these equations some remarks can be made. First, in the special 
configuration M2 ~ ±/i an apparent singularity might arise. Ref. [B] pinpointed this configuration 
which can induce a large f^-scheme dependence in the counterterms 5Mi^2 and 5^, and therefore 
to the annihilation of the LSP into Ws, for a mixed LSP, see also Second, the choice of m^o 
as an input parameter is appropriate only if the lightest neutralino is mostly bino or if the bino 
like neutralino is not too heavy compared to other neutralinos. It is however very easy to switch 
to another scheme or choice of input parameters in the chargino/neutralino sector. For instance if 
the bino like neutralino is the NLSP with mass m^o , like what occurs in the wino scenario that we 
study in this paper, we simply take m^o as input in which case 5Mi is obtained from Eq. Q by 
Jm^o — > 5m^o and Nij — > N2j. 

v) Finally diagonal field renormalisation is fixed by demanding the residue at the pole of the prop- 
agator of all physical particles to be unity, and the non-diagonal part by demanding no-mixing 
between the different particles when on shell. This is implemented in all the sectors. 

vi) Dimensional reduction is used as implemented in the FFL bundle at one-loop through the 
equivalent constrained dimensional renormalisation[27j . 
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2.4 Infrared divergences 

For the processes XX ~^ X,Y — W^^Z", /, .., we can deconipose the one- loop amplitudes in a 
virtual part Mfi^p (for co- annihilation processes with external quarks we also need to add A^^oop) 
and a counter-term contribution A4ct- The sum of these two amplitudes must be ultraviolet finite 
and gauge independent. Due to the virtual exchange of the massless photon and gluon, this sum can 
contain infrared divergencies. This is cured by adding a small mass to the photon and/or gluon, 
and Xg. This is a valid regularisation, even for QCD, for all the processes we study here where the 
non-Abelian nature of QCD does not show up. This mass regulator should exactly cancel against 
the one present in the final state radiation of a photon(gluon). The QED(QCD) contribution is 
therefore split into two parts : a soft one where the photon(gluon) energy E^^g is integrated to less 
than some small cut-off kc and a hard part with E^^g > kc- The former requires a photon(gluon) 
mass regulator. Finally the sum MfZp'^^'^ + ^ct + M'^°l\E^^g < kc) + M^°;^'^{E-f^g > k^) 
should be ultraviolet finite, gauge invariant, not depend on the mass regulator and on the cut k^. 
We take the strong coupling constant at the electroweak scale ag=as[M'^)—Q.l\%. 

2.5 Checking the result 

i) For each process and set of parameters, we first check the ultraviolet finiteness of the results. 
This test applies to the whole set of virtual one-loop diagrams. The ultraviolet finiteness test is 
performed by varying the ultraviolet parameter Cjjv = 1/e, e is the usual regulator in dimensional 
reduction. We vary Cuv by seven orders of magnitude with no change in the result. We content 
ourselves with double precision. 

ii) The test on the infrared finiteness is performed by including both the loop and the soft 
bremsstrahlung contributions and checking that there is no dependence on the fictitious photon 
mass A-y or gluon mass \g. 

iii) Gauge parameter independence of the results is essential. It is performed through the set of the 
eight gauge fixing parameters defined in Eq. (j2.2p . The use of the eight parameters is not redundant 
as often these parameters check complementary sets of diagrams. It is important to note that in 
order to successfully achieve this test one should not include any width in the propagators. However 
we encountered a W boson resonance for the calculation of xX ~^ QQ' ^-nd we had to include a width 
to the W propagator to avoid numerical instabilities; nevertheless this has been done only for the 
evaluation of the hard emission part and not for the virtual and soft part. This will be discussed 
at more length in due course. 

iv) For the bremsstrahlung part, the soft component is added to the virtual corrections and, for 
the hard one, we use VEGAS [28] adaptive Monte Carlo integration package provided in the FFL 
bundle and verify the result of the cross section against CompHepf29]. The hard part is also the 
trickiest, especially when threshold or resonances are encountered as stated above, so for some 
calculations we use BASES [30] provided in the GRACE package[3T] which have a better treatment of 
singularities |20|. We choose kc small enough and check the stability and independence of the result 
with respect to fee- 
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2.6 Effective weighted cross sections 

All cross sections aij where i,j label the annihilating and co-annihilating DM particles i,j can, 
in general, be expanded in terms of the relative velocity Uy , which for neutralino annihilation is 

V — 2/3 = 2./I — 4m?o/s- Away from poles and thresholds, it is a good approximation to write 

y 

CijWy = Gij + bijV^, keeping only the s-wave, Uij, and p-wave, bij coefficients. With T being the 
temperature, x = m^o /T, the thermal average gives 

{(^ij Vij) = aij + 6(6.y - aij/4)/x. (4) 

With (?o = 2 the neutralino LSP spin degree of freedom (sdof), the co-annihilating particle of sdof 
Qi and mass mi contributes an effective relative weight of 

Si.e// = — (1 + Ato^)'^/^ exp(— a; ArTii), Ami — {nii ~ m^o) / rriyO . (5) 
9q ^ ^ 

The total number of sdof is 5e// = ^i9i.eff- An approximation to the relic density is obtained 
through a simple one dimensional integration 

nu2 ( 10 a;;. \ 0.237 X 10-26cm3.s-i ^ p. , ,,3 
\V9*{xf) 24 y XF J J^JF 

Me// - E K^^.)- (6) 



9eff 



aij, bij that are needed to compute <7ij in Eq. ^ are given in cm"^s~^. xp represents the freeze-out 
temperature. g^,{xF) is the effective degrees of freedom at freeze-out. is tabulated inmicrOMEGAs 
and we read it, as well as xp, from there. The weight of a channel (see the percentages we will 
refer to later) corresponds to its relative contribution to J. We find it instructive to consider the 
weighted cross section 

2 (') 

By doing this we somehow normalise the contributions of, in particular, the co-annihilation cross 
sections which can be very large compared to the annihilation cross sections, but which at the end 
do not contribute as much because of the Boltzmann factor, exp(— x Ami). In our plots the weight 
and statistics factors are chosen at freeze-out with x = xp, see Eq. ([5]), and for ease of notation we 
drop the label e//. Since x = Xp is the lowest value of x, see Eq. that enters the calculation 
of the relic density, the weight factor tends to enhance the real contribution of the co-annihilation 
channels. The correct overall weight is in our case given by micrOMEGAs. This fact should be taken 
into account when we compare the figures where the weighted cross sections are shown and the 
tables where the overall weight (extracted from micrOMEGAs) is given. 

Let us stress once more that in order to derive the relic density we do not rely on the approximations 
given in Eqs ([4l[5l[6|) but pass all the cross sections to micrOMEGAs. 

2.7 Interfacing SloopS with micrOMEGAs 

In order to evaluate the relic density, we interfaced SloopS with micrOMEGAs to take full advantage 
of its powerful features concerning the cosmology related part (solving the Boltzmann equations 
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with CO- annihilation, calculation of the effective degrees of freedom, thermal averaging,..). The 
connection between the two codes is summarized in the following: 

1. The MSSM default directory of micrOMEGAs uses SuSpect \3T . In so doing it inherits some 
of the radiative corrections in particular in the spectrum (mass) calculation used in SuSpect . 
From the corrected spectrum micrOMEGAs works out new effective tree-level underlying pa- 
rameters so that gauge invariance is maintained. For the interface we have removed this 
default option of reading from SuSpect and created a model file based on the same tree- level 
lagrangian as the one used in SloopS. In so doing both SloopS and micrOMEGAs calculate the 
same tree-level cross sections. This is also a check on our tree-level cross sections. 

2. The one-loop cross sections of SloopS, appearing into the form of tables showing the cross 
section as a function of the relative momentum p, a{p), are then interpolated and passed to 
micrOMEGAs which substitutes these new corrected cross sections to the corresponding tree- 
level cross sections. Processes that are not corrected (and hence are not substituted) are of 
course kept in the list of processes for the evaluation of the relic density. 

2.8 Finding scenarios in the MSSM parameter space 

The latest limits on the relic density coming from WMAP five years data give the 2a range [1], 

0.098 < n^h^ < 0.122. (8) 

When > mw, channels with gauge bosons in the final state open and LSP's are annihilating 
very efficiently, making it difficult to obtain an absolute value of the relic density within the WMAP 
bounds. In the MSSM this is realized with a neutralino which is mostly wino like or higgsino like 
and its corresponding mass must be around 1 TeV for the latter and 2 TeV in the wino case to be 
in the cosmologically interesting region. Keeping this in mind we did not worry too much about the 
value of the relic cosmic abundance and, instead, we restricted ourselves to get gauge bosons in the 
final state to study the origin of large corrections, if any. Moreover, since the impact of radiative 
corrections can be large, there is not so much sense in picking up a model based on its agreement 
with the current data on the basis of a tree-level calculation and finally we could argue that this 
agreement can be obtained with non-thermal dark matter production, or any other mechanism 
which could avoid too much depletion. This said we also wanted to have a rather light spectrum as 
a supersymmetric solution to the hierarchy problem requires a relatively light LSP and in order to 
have scenarios testable at colliders. Regardless of these remarks we used micrOMEGAs as a guide, 
being careful about translations of effective couplings and input parameters. 

Last but not least, it is important to stress that we did not apply radiative corrections to all sub- 
processes but only to the ones contributing more than 5% to the relic density, the remaining ones 
were included only at tree-level. Most often the processes that we do not correct add up to more 
than 25% of all the processes contributing to the relic density, even if individually their weight is 
small. Therefore when calculating the correction to the relic density, the one-loop corrections we 
compute can get diluted especially if some cancelations occur at one-loop between the processes we 
consider. This point should be kept in mind when we quote the one-loop corrected relic density. 
Ideally we should have corrected all cross sections. This could of course be done with our code 
SloopS and interface to micrOMEGAs , however in these exploratory investigations our aim is to 
uncover the salient features of the radiative corrections to annihilation and co-annihilation of dark 
matter in supersymmetry. 
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For all the scenarios we will give below, the low energy tree-level input parameters are defined at the 
electroweak scale and are: Mi the U{1) gaugino mass, M2 the SU{2) counterpart, ^ the Higgsino 
"mass", M3 the gluino mass, Ma the mass of the pseudoscalar Higgs boson and tan/3. When not 
specified we will take a common sfermion mass. The sfermion trilinear parameter Af is set to zero 
for all generations. We do not impose any gaugino mass unification at the GUT scale. Now let us 
describe the scenarios we study: 

i) Mixed-bino scenario: Usually, assuming gaugino mass universality at the GUT scale leads 

to a bino like LSP. This gives a relic density which overdoses the universe. Relaxing this 
assumption by adding a substantial wino component one can match the WMAP range thanks 
to the opening of gauge boson channels and co-annihilation with . This is easily achieved 
with Ml ~ M2. In our fist scenario Mi ^ 100 GeV so that the mass of the LSP is around 100 
GeV. 

ii) Higgsino scenario : A mainly higgsino neutralino of mass > raw (here ruj^o ~ 230 GeV) 
will automatically annihilate dominantly into gauge bosons and, because of the degeneracy 
with the lightest higgsino like chargino, co-annihilation takes place also. 

iii) Light-wino scenario : A simple way to obtain gauge bosons in the final states of annihilating 
neutralinos is to increase its SU (2) type coupling by decreasing the value of M2 in order to 
have a wino like neutralino whose mass m^o is taken around 200 GeV in this case. Once again 
we have co-annihilation of the Xi with the xt because of the small mass gap between them. 

3 A mixed bino scenario 

The first example we examine corresponds to a neutralino LSP which has a substantial bino com- 
ponent. It is known that an almost pure bino does not annihilate enough to give the right relic 
density in the radiation dominated standard scenario. As it couples mostly to particles with largest 
hypercharge, the R-sleptons, one can increase the LSP annihilation rate by lowering the R-sleptons 
mass. This is typical of the so-called bulk region of mSUGRA. One can also rely on co-annihilation 
with the next-to-lightest supersymmetric particle (NLSP) to reduce the relic abundance. An exam- 
ple is co-annihilation with the f. We have studied these scenarios in [6] including one-loop effects. 
Another solution is to add some wino component by fixing M2 close to Mi, hence the LSP/NLSP 
will have strong couplings with the W boson which will dominate the annihilation rate. This is 
the case we study here. We take chargino masses within the LEP limits while all other particles 
(sfermions, other neutralinos) are heavy. The underlying parameters of the models are given in 
Table □ 



Parameter 


Ml 


M2 


M 


tp Ms 




A, 


Mao 


Value 


110 


134.5 


-245 


10 600 


600 





600 



Table 1: Mixed-bino scenario: Values of the first SUSY set of input parameters. Masses are in 
GeV. 



Mi^2, IJ- are reconstructed from m^o and to-±,to-±. The relevant physical masses are rn^o = 
106.9 GeV, m^± = 124.6 GeV and m~+ =='274.8 GeV. At tree-level m^o = 125.3 GeV. The 
neutralino composition is: Xi 0.945 - 0.20M^ - 0.27H^ - O.IOH^^, where B,W,H^2 denote 
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the U{1) gaugino or bino, the SU(2) gaugino of wino and the higgsmo respectively. The wino 
component is not neghgible. As a consequence annihilation into gauge bosons is dominant. The 
main process XiXi ~* W~^W~ contributes 44% to the relic density. The important co-annihilation 
channels are Xixt Z^W~^ , X1X2 ~^ W^W~ both contributing 5% and XiXi" ~^ ""^ contributing 
8%- XiXi" ~* cs contributes as much as the ud final state. In the following we will refer to only one 
of these quark final states, of course both are counted for the calculation of the relic density. For 
the the XiXi^ co-annihilation, the s-channel exchange of a boson is dominant, see also [33] . 
Before showing our results let us comment on a technicality related to the contribution of the hard 
bremsstrahlung contribution. This concerns the radiative process XiXi" ~^ ud'^, see Fig.[T] 




Figure 1: Real photon emission leading to W return. 



As explained above, when charged/colored particles are involved in the initial/final state, ini- 
tial/final state radiation should be incorporated to have an infrared safe cross section. This emission 
is split into two pieces, soft and hard, and the cross section must be independent of the cut where 
these two parts are defined. Calculating the real emission is a tricky task, especially here. Indeed 
hard photon emission leads to W^-return, bringing the intermediate W on-shell, and therefore would 
lead to numerical instability if no width, Tw^ is provided for the internal W. We have dealt with 
this problem by providing a width to the W only when the radiation is hard, as needed. For the soft 
part no width is introduced in order to achieve the cancelation of the infrared divergence between 
the soft bremsstrahlung and the virtual correction where all masses are real. In summary since 
the W resonance turns up at an energy much larger than the cutoff energy k^, where the matching 
between the soft and hard emission is done, we decided to split the integration on the hard photon 
phase-space in the 2^3 process as follows: 

i) from kc < E^,g < j^{s — (Af^ -I- 2Ty[rMy[r)), no width is implemented 

ii) from 2^(^ ~ i^^w + 2r^M^)) < E^^g < 2^(s — -^^uj) '^i*^ ^ width to the W propagator. 

One must note that for ii) the hard emission is in fact already included in the tree-level process 
XiXi ~^ W^j with the W decaying into a ud{, cs) pair. The process XiXi^ ~^ W^j contributes to 
the relic density but we did not added to our list of cross sections to correct as it contributes less than 
5%. To avoid double counting when calculating the relic density at one-loop, we therefore subtract 
from the list of uncorrected tree-level contributions XiX^ ~* W^'f with the proper branching 
fraction into ud, cs. We will encounter this feature for all other scenarios that lead to such a final 
state and we will treat it in the same way. 
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Another point is related to processes which are initiated through X2 co-annihilation, X1X2 ^ 
W^W~ . We are working with a scheme where the input parameters are m^o , m^± , m-± . Therefore 
m-j^o receives a correction at one-loop. In principle the full one-loop amplitude would write 

■A4i-ioop(™°r"'°°'') = Mi-ioopim^i) + Am^o ^^*'°° (m^o) , (9) 

X2 

where Altrocl™^^) tree-level amplitude, m^o and m^o''~'°°^ = m-j^o -I- Sm^^o is the corrected 

mass. We have neglected the second contribution. This is because the correction to m-j^o is less than 
0.3% for all schemes as shown in Tabic [2] When calculating the relic density we should also in 



Masses [GeV] 








Tree Level 


125.3 


258.1 


270.4 


One Loop - A^-r scheme 


125.13 


258.58 




- MH scheme 


125.31 


258.05 


270.65 


- DR scheme 


125.17 


258.46 


270.47 



Table 2: Mixed-bino scenario: One-loop corrections to the chargino /neutralino masses in GeV in 
the scheme m^am^+ for different tp-schemes: A-tt, DR and MH. 



principle use the corrected physical mass, like for example in the Boltzmann factor, however again 
this is negligible. Results for the weighted cross sections at tree-level and at one-loop are displayed 
in Fig. [51 First of all the QCD and EW corrections to the co-annihilation into light quarks add up 
to about 10% and are practically velocity independent, especially the QCD corrections. The full 
0{a) correction to gauge boson production shows the same feature. The dominant annihilation 
channel into W^W~ gets about -)-7% correction, while X1X2 ~^ W~^W~ is slightly larger with 
9%. XiXi ~^ Z^W'^ is small with about 2%. Had we used a running a at Mz some of the largest 
positive corrections would have been absorbed, however our results show that the full corrections are 
necessary in view of the upcoming precision on the extraction of the relic density. For this scenario 
the corrections, in the a{Q) schemes are positive for all processes we have considered. Nonetheless, 
since the wino component is important in the evaluation of the cross sections because of the SU (2) 
quantum numbers of the final states, we expect these results to be sensitive to the tan /3 scheme, 
since tan (3 enters the mixing of the bino and the wino. 

Table [3] gives in particular the tan/3 scheme dependence. As expected the dependence is not negli- 
gible in particular for the annihilation channel with both LSP in the initial state. The dependence 
weakens for the co-annihilation channels where only one LSP takes part. The Mh scheme is once 
again a bad choice showing once again[6l[T8l[T7] very large corrections. The difference between the 
Arr and DR is about 2% for the co-annihilation channels and 4% for the annihilation channels. 
At the end taking into account the one-loop corrections only to those processes we studied, which 
represent 70% off all processes, the correction on the relic density is about —3% in the A-rr scheme 
and with a defined in the Thomson limit. 
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Figure 2: Mixed-bino scenario. The left/right panel shows the main gauge boson/quark production 
cross sections respectively. All the cross sections are normalised with the corresponding effective 
degrees of freedom given by Eq. ^ with xp — 25.0. Results are shown for the Att scheme o/tan/3. 
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Tree 




DR 


MH 


XiXi - 




[44%] 


a 


+0.81 


+7.6% 


+12.16% 


+29.6% 








b 


+1.219 


+0.78% 


+7.1% 


+24.2% 


yVy^" - 


ud [8%] 




a 


+15.61 


+7.2% 


+9.8% 


+18.8% 








b 


-5.81 


+5.7% 


+8.3% 


+17.4% 


x'Ixt - 


Z"VK+ 


[5%] 


a 


+8.26 


+2.9% 


+4.4% 


+9.7% 








b 


+ 1.42 


-7.3% 


-3.3% 


+10.7% 






[5%] 


a 


+17.81 


+9.0% 


+11.1% 


+18.2% 








b 


+11.86 


+4.8% 


+7.3% 


+16.1% 










0.108 


0.105 


0.102 


0.097 












-2.8% 


-5.6% 


-10.2% 



Table 3: Mixed-bino scenario: Tree-level values of the s-wave (a) and p-wave (b) coefficients in 
units 10~^^cm^s~^, as well as the relative one-loop corrections in the A-^r, DR, and MH scheme. 
The percentages in the first column refer to the percentage weight, at tree-level, of that particular 
channel to the relic density. 

4 A light Higssino scenario 

A pure Higgsino state could give an interesting relic density and, as the XiX"^" and Xixf ^""^ are 
large, annihilates mainly into WW and ZZ final states. Besides, as there are three Higgsino like 
states (two neutralinos and one chargino) whose mass splitting is small especially if gaugino masses 
are large , m^o ~ m^o ~ ~ I'^l' co-annihilation between the LSP and the other higgsino states 
is important. With such efficient annihilations the relic density would be small if the Higgsino like 
LSP is too light. Nonetheless it gives favourable prospects for dark matter direct detection. The 
scenario we have chosen is described in terms of the underlying parameters given in Table [H 



Parameter 


Ml 


M2 


fi tp 


Ms 






Mao 


Value 


400 


350 


-250 4 


1000 


650 





800 



Table 4: Higgsino scenario: Parameters defining our higgsino model with little mixing. Masses are 
in GeV. 

The LSP neutralino with mass m^o = 234 GeV has a composition Xi = O.ll-B - 0.31M^ - 0.70.^^* - 

0.63-ff2 1 indicating it is dominantly a higgsino state. Co-annihilation between the Xi and occurs 
since m-± — 242.9 GeV. All other particles, including y^, are heavy enough and therefore do not 
take part in the co- annihilation. As the higgsino component of the neutralino is important, it will 
couple mostly to the W and, compared to the mixed-bino case, to the boson also through the 
XiX2^° coupling. 

Dominant tree-level processes relevant for the computation of the relic density are the same as in 
the previous mixed-bino case except for the co-annihilation between the first two neutralinos which 
is Boltzmann suppressed due to their larger mass splitting and smaller couplings. The dominant 
processes are XiXi ~^ W^W~ contributing (at tree-level) 26% to the relic density, XiXi^ ~^ ud{cs) 
with 12%(12%), XiX? Z°Z° with 9% and XiX^" ^ Z°W+ with 6%. 
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Our results for the cross seetions both at tree-level and at one-loop are displayed in Fig. [31 They 
are shown for A^-r scheme of tan (3 and where the input for a is in the Thomson limit. 
Compared to the mixed bino case, the QCD and EW corrections to co-annihilation into light quarks 
are smaller and no cancelation between the two occurs. The overall correction is almost velocity 
independent and ranges between 10% to 8%. The corrections to gauge boson production are smaller 
for x?X? ^ W+W- and XiXt ^ Z°W+ and about 10% for XiXi ^ Fig. [3] shows a very 

interesting dynamical effect in the one-loop correction to XiXi ~^ W^W^ , for ~ 0.3. The 

bumps are in fact due to the opening of the threshold in the loop, as can be checked explicitly 
for this value of the velocity and the mass of the LSP compared to that of the . This dynamical 
structure can not be described by a simple a + bv^ parametrisation of the cross section. Compared 
to the bino case we have studied in the previous section the tan/3 scheme dependence is small. 
The dependence is shown in Table [5] where we also give the results in terms of the s-wave and 
p-wave coefficients for a fit in the range v'^ < 0.3 so that we avoid the dynamical structure we 
have just pointed at. The difference between the Arr scheme and the DR scheme is below 1% for 
all processes we studied, while the MH scheme gives larger corrections but within 2% compared 
to the DR. For quark production the scheme dependence is even negligible. The overall 0{a) 
corrections in this scenario, though not negligible, are not that large with a defined in the on-shell 
scheme in the Thomson limit. Moreover corrections coming from boxes and final state radiation 
are often dominant. This suggests that to grab most of the radiative corrections a full calculation 
is needed. Within our approach of not correcting the processes that contribute less than 5% to the 
relic density, the processes we have considered contribute in total only 65%. In this approach we 
find a correction to the relic density of —2.5% in the Att scheme and —2.4% in DR. In the MH 
scheme the correction is little changed to —3.3%. 









Tree 




DR 


MH 




^ W+W- [26%] 


a 


+ 11.84 


+4.3% 


+5.1% 


+6.8% 






b 


-f4.17 


+ 12.7% 


+13.4% 


+14.9% 


XiX^ - 


^ ud [12%] 


a 


+ 15.28 


+6.8% 


+7.0% 


+7.3% 






b 


-5.31 


+30.4% 


+30.7% 


+31.3% 


x^x'i' - 


^ [9%] 


a 


+4.28 


+ 10.4% 


+9.6% 


+7.8% 






b 


+1.83 


+ 12.7% 


+12.0% 


+10.5% 


XiX^ - 


Z"W+ [6%] 


a 


+6.99 


+1.7% 


+2.1% 


+2.9% 






b 


-0.51 


+85.6% 


+86.5% 


+88.4% 








0.00931 


0.00909 


0.00908 


0.00904 


sn^h^ 








-2.4% 


-2.5% 


-3.3% 



Table 5: Higgsino scenario: Tree-level values of the s-wave (a) and p-wave (b) coefficients in units 
lO^^^cm'^s"^ in the higgsino scenario, as well as the relative one-loop corrections in the Arr, DR, 
MH scheme. The percentages in the first column next to the process refer to the percentage weight, 
at tree-level, of that particular channel to the relic density. The fit into a and b is done in the range 
<v^ < 0.3. 
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Figure 3: Higgsino scenario: The left/right panel shows the main gauge boson/quark production 
cross sections respectively. All the cross sections are normalised with the corresponding effective 
degrees of freedom given by Eq. ^ with xp = 27.6. 
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5 A light wino scenario 



Models with a light wino as the dark matter candidate occur in theories like AMSB[34], string 
compactifications(3S] and also split-supersymmetry[351[37]. The advantage of a light wino is that it 
has a large annihilation cross section, relevant for indirect detection, but the main drawback is that 
it predicts a small thermal relic abundance in the standard cosmological scenario and non-thermal 
production has to be invoked to recover the correct relic density. The underlying parameters of the 
model are given in Table [S) 



Parameter 


Ml 


M2 






M3 








A, 




Value 


550 


210 


-600 


30 


1200 


387 


360 


800 





700 



Table 6: Light-wino scenario: Values of the fourth SUSY set of input parameters. Masses are in 
GeV. 

The LSP is now essentially wino with a composition Xi = 0.005S - 0.99W - 0.15H^ ~ O.Q5H^ 
and mass m^o = 206.6 GeV. The LSP is highly degenerate with the Xi , their mass difference 
is Am ~ 0.05 GeV and consequently sizeable co-annihilations occur in the determination of the 
relic density. With so small mass difference, co-annihilation processes are important. Products of 
annihilation/co-annihilation processes are into gauge bosons (and some light quarks). The domi- 
nant processes are the following: XiXi ^ W+W- [13%], xtxt_ ^ W+W+ [12%], Xixt ^ Z"W+ 
[12%], xtxl Z°Z" [7%], xtxl W+W- [7%], x°ixt ^ [7%]. Note in passing that we have 
taken a large value of tan (3. 

Before we present our results, a word about the renormalisation scheme and the choice of input 
parameters especially as concerns the neutralino/chargino sector is in order. The results we will 
show are based on taking the mass of the LSP as input (beside the masses of the charginos) . One 
might argue that this is not optimal in order to reconstruct the system Mi , AI2 , /i, especially for 
extracting Mi which is sensitive to the bino-component. One might even expect that at one-loop this 
scheme would not be suitable since the A^n element of the orthogonal matrix in the neutralino sector 
is very small leading to a large contribution from the counterterm SMi , see Eq. ([3]) . For this reason 
we have been careful in also taking the scheme where the input parameters are (ruj^o , rn^+ , w -+ ) . 

Fortunately, as we can see in Fig. |4]for the process XiXi ~^ W+W~ , the difference between the 
two schemes is within less than 0.4%. We have checked that for other processes in this scenario the 
difference is also negligible. 

The tree-level cross sections and the full one-loop corrections are shown in Fig. [5] as a function 
of the relative velocity. We note that, at tree-level, these cross sections are s-wave dominated. 
For bosonic final states the velocity dependence is modest, compared to the co-annihilation into 
quarks. For the latter, the electroweak and QCD corrections to Xixt ud are relatively large, of 
order 30%, but they almost cancel each other. The overall correction is about -1-5% and practically 
independent of the velocity. The annihilation process and the Xixf co-annihilation processes show 
an almost constant correction of order 10%. The co-annihilation processes show an interesting 
behaviour in the case where both co-annihilating particles are charged, the cross sections reveal a 
very large correction at very small relative velocity. This correction is the one-loop manifestation 
of the non-relativistic Coulomb- Sommerfeld effect [19] . With the tree-level cross section denoted 
as (To E^nd uqv = + b^v^ , at vanishing relative velocity the one loop cross section for chargino 
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Figure 4: Relative difference for XiXi W~^W~ between the scheme with (m^o , m -+ , to -+ ) and 
(m-^o, m~+, TO^+) as function of relative velocity, tp is within the A^-t scheme. 



annihilation, cr^^°°^ is such that 

l-loop f iTTQ r ~±~T 

" I for xf xf 

We thus expect the one-loop cross section cti to be 

aiv = ax + biv'^ + iraci/v with ci = zbao- (10) 

To exactly quantify the Sommerfeld effect in our calculation, we have also fitted the one-loop cross 
section in the form of Eq. p^ . An example of such an exercise is given in Fig. [6] for Xi^X^ ~^ 
W~^W^. We see therefore that our calculation captures this effect extremely well, indeed we 
obtain here that ci/ao 1.015 which indeed very close to the analytical result, ci/ao = 1. This 
is important because this effect needs to be summed up to all orders. In our approach we will 
therefore subtract it from the one-loop correction and replace it by the resummed all order result in 
the final result. The result of this subtraction on the processes xtxi ^ Z^Z^, Xi^Xi^ ~* W^W^ 
and x^xr ~^ W^W~ is shown in Fig. [7] 

As one can see once the Coulomb-Sommerfeld contribution is removed, one is left with a smooth 
cross section which is almost velocity independent. 

Looking carefully at the results for Xi'xT ~^ W^W~ we note that there is still a slight increase at 
small V. This is a residual effect of the weak Sommerfeld contribution, see Fig. [H mediated by a 
charged W that it is noticeable even for a not too heavy chargino. In fact a similar effect is also 
present in XiXi ~^ W^W~ . Zooming in on the region of small relative velocity we see a kink, see 
Fig. [51 around y/s ~ 413.3 GeV corresponding to u ~ 0.04 which corresponds to the opening of the 
threshold of chargino production. In Xi^xT ~* W^W~ we only see the tail of the opening of the 
threshold. 

Apart from these interesting but numerically small features, let us mention that the tp scheme 
dependence is negligible, it is below 0.1%. Our results show that corrections to the individual cross 
sections can get large at all relative velocities, even after subtracting the QED Sommerfeld effect. 
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Figure 5: Light-wino scenario: The left/right panel shows the main gauge boson/quark production 
cross sections respectively. All the cross sections are normalised with the corresponding effective 
degrees of freedom given by Eq. ^ with xp = 29.9. 
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Figure 6: Lightwino scenario. Fits to the s-wave, p-wave and Sommerfeld factors for Xi Xi 
W^W^ . In particular note that the fit in the parametrisation of Eq. ilO\) gives Ci/ao = 1.015. 
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Figure 7: Light-wino scenario: Results for one-loop corrections in the A^-t scheme where the QED 
Sommerfeld effect has been subtracted. 



For example, Xi^xT ~^ W^W~ gets about +30% correction, while both XiXi —> W^W^ and 
XiXi" ~^ Z'^W'^ get a —10% correction. The dominant cross section XiXi ~^ W^W^ receives 
a correction of about +15%. The corrections for the other processes are more modest. These 
corrections are within the on-shell scheme by using q;(0). We see that had we used a{M'^) the 
correction to XiXi ~^ W'^W~~ would be small, but this is not true for most of the other cross 
sections where genuine corrections, including hard radiation effects are important and must be taken 
into account. This said, when we combine all the cross sections, taking into account their statistical 
weight, substantial cancelations occur between the different contributions. Add to this that the 
cross sections we have considered account for about only 65% of all cross sections contributing to 
the relic density, since we have not considered those contributing individually less than 5%, we find 
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a quite modest (within the a(0) scheme) correction to the rehc density of about —2%, see Table [71 
This full one- loop correction is practically unchanged if we instead resum the 1/v Sommerfeld 
effect. A similar result was found when we studied f co- annihilation [B]. This is due to the fact that 
temperature effects provide a cut-off and the l/v enhancement is tamed after thermal averaging, 
(X /o°°(dz; v"^ e~" /^^) {crv). Our results are summarised in Table [T] The results are presented in 
terms of the s-wave and p-wave coefficients as well as the Sommerfeld l/v coefficient. The correction 
to the relic density is given for the full one-loop, including the one-loop l/v threshold correction, as 
well as after resumming the l/v contributions. Another word of warning about the interpretation 
of the corrections in terms of the s-wave and p-wave coefficients (a and b). The corrections to the 
p-wave coefficients may seem very large here, especially if the corresponding correction to a is large. 
This is not an indication that the radiative correction on the total cross is very large. Indeed all the 
cross sections here are s-wave dominated, so that the correction on the s-wave is a good measure 
of the total correction and when modulated with the statistical weight gives a good approximation 
to the correction on the relic density. 
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+46.49 
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-104.4% 
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+53.34 








0.00215 


0.00211 




1-loop 






-1.9% 




with resum. Sommerfcld 






-1.9% 



Table 7: Light-wino scenario. The table summarises the effect of the full order corrections on the 
dominant processes that contribute more than 5% to the relic density. The relative contribution is 
given in [ ] next to the process. The tree-level cross sections are given through the fit av — a + bv"^ 
in the range < f ^ < 0.3. At the one-loop level, The fitting function is then av = a-\- bv"^ + CT:a/v. 
The coefficients a, b and c are given units o/ 10~^^cin'^s~^. The relic density is calculated by taking 
into account all other processes, which however are not corrected at one-loop. The Table also gives 
the correction after summing the 1/v QED contribution at all orders. As the tp-scheme dependence 
is less than 0.1%, only one tfi scheme A^-t is presented. 



6 Conclusions 

Extraction of the relic density will soon provide a measurement of this quantity at the 1% level. On 
the theoretical side one must therefore provide a prediction which is at least as precise. In particular, 
if the particle physics component in terms of computation of the annihilations and co-annihilations 
cross sections are under control, one can indirectly test the cosmology of the Universe. With this 
precise measurement we can even gain insight into the particle physics model that could be combined 
with measurements at the colliders. The work in this paper continues the investigations we have 
made [6.^ concerning the impact of the radiative corrections on the annihilation and co-annihilation 
cross sections of a neutralino dark matter in the MSSM. Here the emphasis is on processes with 
dominant gauge boson production channels. We have considered three models with relatively light 
LSP in the 100-200 GeV range, i) a dominantly bino with some admixture of wino, ii) a higgsino 
like and a iii) wino like LSP. Our study shows that it is not easy to find a general common feature 
of the radiative corrections. For example, within the same tan (3 scheme and for relative velocities 
relevant for the evaluation of the relic density, the dominant process XiX? W^W~ gets about 
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15% correction in the wino case but only 4% in the higgsino case, while the bino is in between. 

Also the corrections wc have just quoted arc within a scheme where the electromagnetic constant 
is defined at low energy in the Thomson limit. A naive use of a(M|) would suggest that most 
of the corrections in the dominant process in the wino case is absorbed, but this would not be 
true for the other processes nor for the the other scenarios. This still docs not take into account 
the effects of final state radiation. For example in the same wino scenario, the 0{a) correction 
to XiXi ~^ W~^W~ is large and reaches about 30%. In general the corrections to the different 
contributing processes for the same scenario can be quite different, in the case of the wino the 
overall effect on the relic is a cancelation of the corrections between the different processes. With 
this in mind and the fact that we did not correct processes that, individually, contribute less than 
5% to the relic, we find that the overall effect on the relic is small, —2% in the on-shell scheme 
with a in the Thomson limit. In this paper we have also pursued our investigation of the effect of 
the tan/3 scheme dependence on many observables, not necessarily dark matter annihilation. We 
confirm once more that the M H scheme is not appropriate while DR and A^r give generally similar 
results. In the scenarios we have studied, in fact the scheme dependence is an issue only for the bino 
case. This could have been expected as the bino couples to W^s only through mixing where tan/3 
is important. We have also uncovered in the case of co-annihilation electromagnetic Sommerfeld 
effects for vanishingly small relative velocity. However the result of the full one-loop and that of 
resumming this effect is not noticeable at the level of the relic density evaluation, thermal averaging 
smoothes out the effect. Althoiigh the effect is numerically quite small, in the case of the wino we 
noticed the effect of the electroweak Sommerfeld enhancement. This will become more important 
for higher wino masses that we will study in another paper. To sum up, it is important to stress 
that we now have the tools to perform automated calculations relevant for a precise evaluation of 
the dark matter annihilation cross sections. 
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